Non-Fermi-liquid phases in the two-band Hubbard modeh 
Finite-temperature exact diagonahzation study of Hund's rule coupling 
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The two-band Hubbard model involving subbands of different widths is investigated via finite- 
temperature exact diagonahzation (ED) and dynamical mean field theory (DMFT). In contrast 
to the quantum Monte Carlo (QMC) method which at low temperatures includes only Ising-like 
exchange interactions to avoid sign problems, ED permits a treatment of Hund's exchange and 
other onsite Coulomb interactions on the same footing. The role of finite-size effects caused 
by the limited number of bath levels in this scheme is studied by analyzing the low-frequency 
behavior of the subband self-energies as a function of temperature, and by comparing with 
numerical renormalization group (NRG) results for an effective one-band model. For half-filled, 
non-hybridizing bands, the metallic and insulating phases are separated by an intermediate mixed 
phase with an insulating narrow and a bad-metallic wide subband. The wide band in this phase 
exhibits different degrees of non-Fermi-liquid behavior, depending on the treatment of exchange 
interactions. Whereas for complete Hund's coupling, infinite lifetime is found at the Fermi level, 
in the absence of spin-fiip and pair-exchange, this lifetime becomes finite. Excellent agreement 
is obtained both with new NRG and previous QMC/DMFT calculations. These results suggest 
that-finite temperature ED/DMFT might be a useful scheme for realistic multi-band materials. 

PACS. 71. 20. Be Transition metals and alloys - 71.27-fa Strongly correlated electron systems 



1 INTRODUCTION 

Strongly correlated materials exhibit a wealth of fas- 
cinating physical phenomena associated with complex 
single-electron and many-electron interactions. Transi- 
tion metal oxides, for example, tend to have partially 
filled shells of highly correlated d electrons, surrounded 
by complicated lattice geometries, with many atoms and 
electrons per unit cell. The theoretical description of 
the electronic properties of these materials is a challeng- 
ing topic in condensed matter physics. Significant ad- 
vances were achieved during recent years via the dynam- 
ical mean field theory (DMFT) 0, which provides a 
treatment of single-electron and many-electron interac- 
tions on the same footing. 

For realistic materials, DMFT has been used exten- 
sively in combination with the quantum Monte Carlo 
(QMC) method 0. This approach has the advantage 
that it can be readily applied to systems consisting of two 
or more subbands. It has the drawback, however, that, 
to avoid sign problems at low temperatures, it includes 
only Ising-like exchange interactions 0. Improvements 
of the QMC method to include the full Hund's coupling 
are so far limited to T = or rather high temperatures 
|M IH Q' In view of this limitation there is clearly a 
need to explore alternative methods that are applicable 
to multi-band materials and complete onsite exchange 
interactions. 

The aim of the present work is two-fold: First, ex- 
act diagonahzation (ED) 0, is proposed as a poten- 
tially highly useful impurity solver for finite-T DMFT 
studies of realistic systems. The attractive feature of 
this approach is that, in contrast to finite-T QMC, on- 



site Coulomb and exchange interactions are treated on 
the same basis. Second, we apply finite-T ED/DMFT 
to a highly nontrivial system which has recently re- 
ceived considerable attention, namely, the two-band 
Hubbard model consis ting of subbands of different widths 

iiiEiiiiiiiiiiiliaiiiiiiiiiiiiiiaia. The 

phase diagram of this system was recently evaluated in 
Refs. Igj ,17.] . Here we focus on the electronic proper- 
ties of the so-called orbital-selective phase in which the 
narrow band is insulating while the wide band is still 
metallic. To examine the infiuence of finite-size effects in 
the ED approach, and to investigate the electronic prop- 
erties in the limit of low temperatures, we have also per- 
formed numerical renormalization group (NRG) DMFT 
calculations for a simplified two-band model which is par- 
ticularly suited for the intermediate phase. 

The main result of this paper is that the two-band 
ED/DMFT calculations provide a correct picture of 
the electronic properties of the Hubbard model involv- 
ing nonequivalent subbands, including the unusual non- 
Fermi-liquid properties of the orbital-selective phase. 
This outcome is remarkable since, for computational rea- 
sons, the number of bath levels per impurity orbital is 
necessarily smaller than in analogous one-band models. 
Nevertheless, despite this limitation, the ED results are 
in good qualitative or, in some cases, quantitative agree- 
ment with the NRG results. In the Ising case they 
also agree very well with previous QMC/DMFT results 

[HIIil 

The overall consistency with the NRG and QMC cal- 
culations suggests that finite-T ED/DMFT might indeed 
be very useful in the future to analyze realistic materials, 
especially, if full diagonahzation is replaced by finite-T 
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Lanczos methods [23, order to be able to deal with 

larger cluster sizes. 

Hubbard models involving nonequivalent subbands 
are relevant for compounds such as Ca2-a;Sra;Ru04 
and Naa;Co02, where, as a result of the nearly two- 
dimensional lattice structure, coexisting narrow and wide 
bands arise naturally. Thus, onsite Coulomb energies can 
be simultaneously large and small relative to the widths 
of important subbands. As a function of doping concen- 
tration, both materials give rise to a remarkably rich se- 
quence of phases, inclu ding superconductivity and Mott 
insulating behavior j^iEHl- Evidently, the competition 
between multiple kinetic energy scales and local Coulomb 
and exchange energies is an important feature of these 
strongly correlated systems. 

A consistent treatment of Coulomb and exchange in- 
teractions in materials of this kind is important since it 
has recently become clear that the Hund's coupling has 
a decisive influence on the nature of the Mott transition 
H H E H [H 113. In fact, the different treatments 
of exchange terms in earlier finite-T QMC 0, llll 
and zero-T ED calculations have given rise to some 
confusion, with apparently contradictory results. As was 
clarified in Ref. [3 for T = and in Ref. [13 for T > 0, 
however, the QMC and ED results are in agreement pro- 
vided that exchange interactions are treated in the same 
manner. Thus, for full Hund's coupling the two-band 
Hubbard model exhibits successive first-order transitions. 
In striking contrast, in the absence of spin- flip and pair- 
exchange only the lower transition remains first-order 

Moreover, as will be discussed in detail below, the na- 
ture of the intermediate phase depends in a subtle man- 
ner on the treatment of exchange interactions. In this 
regard the ED and NRG results yield the following pic- 
ture: For full Hund's coupling, the wide band has infinite 
lifetime at Ep but does not satisfy Fermi-liquid criteria 
at finite frequencies. This finding is consistent with re- 
cent results obtained by Biermann et al. in T = 
two-band ED/DMFT calculations. For Ising exchange, 
instead, correlations are significantly enhanced and the 
lifetime becomes finite even at Ep, in agreement with 
QMC calculations IHEJ. 

Thus, the Mott transition in the Hubbard model in- 
volving different subbands does not consist of equivalent 
sequential transitions. Instead, when the narrow band 
becomes insulating, the wide band is forced into a bad- 
metallic state whose deviations from Fermi-liquid behav- 
ior depend in a qualitative manner on the treatment of 
exchange interactions. 

In the past, finite-T ED/DMFT methods have been 
applied mainly to single-band cases, where the incorpo- 
ration of an appropriate cluster of bath levels is feasible. 
Throughout this paper we consider two impurity levels, 
each surrounded by either 2 or 3 bath levels, giving total 
cluster sizes ris = 6 or ris = 8 per spin. Below we demon- 



strate that even a cluster size of Hs = 6 provides qual- 
itatively correct subband self-energies. For instance, at 
T = 10 meV the critical Coulomb energies of the orbital- 
selective Mott transitions for = 6 differ by only about 
0.1. . . 0.2 eV from those derived for = 8. This finding 
is interesting since it suggests that finitc-T ED/DMFT 
calculations for more realistic three-band models using a 
cluster size of rig = 9 (three impurity levels, each cou- 
pled to two bath levels) should be useful. This would 
allow one to re-examine the Mott transition in systems 
that have been studied previously using QMC/DMFT for 
Ising exchange interactions. 

A more accurate representation of low-frequency prop- 
erties at low temperatures requires three bath levels per 
impurity orbital: one near the Fermi level to provide ad- 
equate metallicity, and two for the upper and lower Hub- 
bard bands, giving = 8 per spin. This extension leads 
to a significant reduction of finite-size effects. Several 
comparisons of results for Us — 6 and Ug ~ 8 are pro- 
vided below to illustrate the range of applicability of the 
two-band ED/DMFT approach. 

The outline of the paper is as follows: In Section 2 the 
multiband Hubbard model is specified and its numerical 
solution via the finite-T exact diagonalization method 
is discussed. Section 3 describes the effective one-band 
model used in the NRG approach to evaluate the elec- 
tronic properties of the wide band in the intermediate 
phase when the narrow band is insulating. In Section 4 
the ED/DMFT is applied to the purely metallic phase 
just below the first Mott transition. Section 5 deals with 
the intermediate phase in the presence of full Hund's cou- 
pling, where wide band exhibits infinite lifetime at Ep, 
combined with non-Fermi-liquid behavior at finite fre- 
quencies. In Section 6 this intermediate phase is consid- 
ered in the absence of spin-flip and pair-exchange terms, 
giving rise to finite lifetime even at Ep. In Section 7 
the two-band ED approach is applied to the case studied 
previously within QMC/DMFT JJj. Section 8 presents 
analogous results for the case considered recently within 
QMC/DMFT by Biermann etaL Section 9 contains 
a brief discussion of iterated perturbation theory (IPT) 
DMFT results for the two-band model and Section 10 
provides a summary and outlook. 

2 MULTIBAND EXACT DIAGONALIZATION 
METHOD 

The two-band Hubbard model for non-equivalent sub- 
bands is represented by the Hamiltonian: 

H = Hq + Hi + H2 

— ^ ^ tl^-jiiCj^^^Cjjiifj 
ImicT 

li laa' 
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H2 ^ -j'^[c+^QiiC+^Q2T +H.C.] 

I 

-•^'Xl[^iiT'^niC;2TC/2i +H.C.] , (1) 



where c^^ and cua are creation and annihilation opera- 
tors for electrons at site I in orbital i = 1,2 with spin 
a and nucr = cf^^cua. H.c. denotes Hermitian conjugate 
terms. is the single-particle Hamiltonian which we 
represent by half-filled, non-hybridizing bands with den- 
sity of states Ni{uj) = 2/{'KDi)y/l — {uo/DiY and widths 
Wi ^ 2 eV, = 4 eV, where = 2D,. Hi repre- 
sents the anisotropic, Ising-like onsite Coulomb and ex- 
change interactions and H2 additional spin-flip and pair- 
exchange contributions. For isotropic Hund's coupling 
one has J' = J = {U — U')/2. Below, we discuss results 
for J' = J ds, well as J' — 0. The latter case corre- 
sponds to the Ising-like exchange treatment in previous 
QMC/DMFT calculations 

The ED/DMFT results are derived from a two-band 
generalization of the approach employed for single bands 
PJ, |3| . The effective Anderson impurity Hamiltonian in- 
cludes impurity levels ei_2 and bath levels ek=^...n^ (see 
Fig. 1). For particle-hole symmetry £1^2 = 0. Each impu- 
rity level interacts with its own bath, so that for ris = 6 
impurity level 1 couples to bath levels 3,4 with £4 ~ —£3, 
while level 2 couples to bath levels 5,6 with £5 — —£5. 
For rt<; = 8 impurity levels 1,2 couple in addition to the 
bath levels £7 = and £3 = 0, respectively. 

Since at T > all states of the impurity Hamiltonian 
are used in the construction of the subband Green's func- 
tions, for ris = 6 the largest matrix to be diagonalized 
corresponds to the sector n-\ = ni — 3, with dimension 
[ns\/{{ns/2)\ff = 400. For 8 this dimension in- 

creases to 4900. We neglect hybridization between bands, 
so that self-energies and Green's functions are diagonal 
in orbital space. Denoting eigenvalues and eigenvectors 
of the impurity Hamiltonian by and the subband 
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FIG. 1: Level diagram for ED scheme. Blue lines: impurity 
levels; red solid (plus dashed) lines: bath levels for ris = 6 
(us = 8). Green lines: Coulomb and exchange interactions 
between impurity levels; red dotted lines: hopping interac- 
tions between impurity and bath levels. 
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FIG. 2: Zi{U) as a function of [/ for J = J' = (7/4 at T = 
10 meV, calculated within ED/DMFT for different cluster 
sizes: (a) Us — 8; (b) ris = 6. Solid (red) curves: narrow 
band, dashed (blue) curves: wide band. 



Green's functions are evaluated from the expression 



E,, - E„ - I 



(2) 



Here, j3 ~ l/fc^T, ojn — {2n + l)7r//3 are Matsubara fre- 
quencies, and Z = exp(— is the partition func- 
tion. I = denotes the impurity site. Since we consider 
only paramagnetic phases the spin index of Green's func- 
tions and self-energies is omitted for convenience. Be- 
cause of the particle-hole symmetry of the present sys- 
tem, Green's functions and self-energies are purely imag- 
inary quantities. 

To provide an overview of the Mott transitions in 
this model and to illustrate the sensitivity of the crit- 
ical Coulomb energies to the cluster size used in the 
ED/DMFT, we show in Fig. El the quantities Z^{U) = 
1/[1 — lm'Si{iuJo) / loq] for Ug ~ 6 and Us — 8 at T — 
10 meV, assuming J' = J — U/A. Here I]i(iwo) are 
the subband self-energies at the first Matsubara fre- 
quency. Thus, for isotropic Hund's exchange and Us — 8, 
Uci ~ 2.1 eV and Uc2 ~ 3.1 eV. (In our notation the 
Uci=i.2 refer here to the subband critical Coulomb ener- 
gies for increasing U and should not be confused with the 
stability boundaries for increasing vs decreasing U at the 
individual Mott transitions.) Near both critical Coulomb 
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energies Zi{U) exhibit hysteresis behavior characteristic 
of first-order transitions. For rig = 6, the common metal- 
Hc phase is stable only up to about Ud ~ 2 eV. Thus, 
inclusion of zero-energy bath levels supports the metallic 
character of the DMFT solution. A shift of about 0.2 eV 
is found for Uc2 where the wide band becomes insulating. 
These shifts are consistent with single-band ED/DMFT 
results [l^ upon increasing Ug from 3 to 4; only minor 
additional shifts occur in this case between = 4 and 
the fully converged results for = 6. 

According to the QMC/DMFT results discussed in 
Ref. the Mott transitions for J = 11/4 and J' = 0, 
i.e., for Ising-like exchange, are located at Ud ~ 2.1 eV 
and Uc2 ~ 2.7 eV. Moreover, apart from the shift of 
Uc2, this transition is no lo nger first-order. The two-band 
ED/DMFT results in Ref. [l^ for = 6 confirmed this 
fundamental difference between the J' = J and J' = 
treatments. 

As will be discussed in detail below, the wide band in 
the intermediate phase, i.e., for Ud < U < Uc2, does not 
satisfy Fermi-liquid criteria. Thus, the fact that Z2{U) 
as defined above is finite in this region does not imply 
existence of ordinary quasiparticles. In the case of full 
Hund's coupling, the imaginary part of the self-energy of 
the wide band vanishes in the low- frequency limit, but 
does not increase linearly at small iwn, implying non- 
quadratic variation at real frequencies. For Ising-like 
exchange, Iml]2(*'^n) remains finite for uj„ 0, im- 
plying finite lifetime even for states at Ep- Plots like 
those in Fig. 1 are nevertheless useful since they permit 
a convenient identification of phase transitions. Thus, 
Zi{U) — >• at Uci indicates that the narrow or wide 
subbands become insulating for U > Ud, respectively. 
(Of course, due to the discrete representation along the 
imaginary frequency axis at finite T, Zi(U) does not fully 
vanish.) 

To check the convergence of the ED results with clus- 
ter size Us we compare them with NRG calculations for 
an effective one-band model which is most suitable for 
the intermediate phase and which is specified in the fol- 
lowing section. Since the narrow band is insulating in 
this phase one of the key features of the effective model 
is the omission of one-electron hopping in the narrow 
band. Accordingly, the upper and lower Hubbard peaks 
will be centered at about ±C//2, but the infiuence of their 
width is neglected. (In a local description, the width of 
the Hubbard bands is approximately given by [ij.) In 
Fig. 131 we show that this approximation indeed has only 
a minor effect on the electronic properties of the metal- 
lic wide band. Keeping W2 = 4 eV fixed and chang- 
ing Wi between 0.1 W2 and 0.5 W2, we notice that, both 
for isotropic and anisotropic exchange coupling, the self- 
energy of the wide band atU — 0.6 W2 = 2.4 eV is almost 
unaffected by the original single-particle width of the in- 
sulating narrow band. The detailed electronic properties 
of these phases will be discussed in Sections 5 and 6. 
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FIG. 3: Imaginary part of self-energy of wide band in in- 
termediate phase for decreasing widths Wi of narrow band, 
with fixed W2 ^ 4 eV, U ^ 2.4 eV, = 8, T = 10 meV. 
Alternating solid (red) and dashed (blue) curves: Wi = 
0.4, 0.8, 1.2, 1.6, 2.0 eV (from bottom), (a) J = J' ^ [7/4; 
(b) J = C//4, J' = 0. 



3 NUMERICAL RENORMALIZATION GROUP 
APPROACH FOR EFFECTIVE ONE-BAND 
MODEL 

In the limit of an insulating narrow band interact- 
ing with a metallic wide band, the two-band model can 
be simplified to an effective one-band model by elimi- 
nating the high energy states associated with the upper 
and lower Hubbard bands of the narrow band 21]. Ne- 
glecting the one-electron hopping in the narrow band, 
i.e., assuming Wi = 0, and fixing the occupations of 
these orbitals, the terms in Eq. (1) involving interorbital 
Coulomb and pair-exchange interactions disappear. The 
effective Hamiltonian then reduces to 

Ima I 

- Y,[^JStst + 2J\S+s- + Sis+)] . (3) 
; 

In this model, the low energy degrees of freedom of the 
narrow band are represented by local moments Si . These 
interact ferromagnetically with the local spin density Si 
of the wide-band conduction electrons which are also sub- 
ject to a local Coulomb repulsion U. The model corre- 
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spends to a double-exchange model with Coulomb inter- 
actions amongst the itinerant electrons or, equivalently, 
to the ferromagnetic Kondo lattice model with interac- 
tions in the band. The antiferromagnetic case, J < 0, 
for U = 0, has already been investigated in the context 
of heavy fermions [23 ■ We adapt these calculations to 
the case of interest here, namely, ferromagnetic exchange, 
J > 0, and U > 0. The equivalence between the full two- 
band model and the effective model holds provided U 
is large enough so that a description of the low energy 
degrees of freedom of the narrow band in terms of local 
moments is possible and provided that J' ^ U/2. The 
latter condition corresponds to eliminated excited states 
of the full model being far from the first excited state of 
the effective model. Unless J is chosen to be very small, 
we see that for realistic values of J « U/A, the equiva- 
lence will hold well for Ising exchange coupling (J' — 0) 
but less well for isotropic Hund's coupling (J' = J). In 
either case, the effective model should be increasingly ac- 
curate in the limit ix), T — > 0. 

Within DMFT, wc need to solve an effective quantum 
impurity model corresponding to a S" = 1/2 Kondo impu- 
rity coupled ferromagnetically with conduction electrons 
subject to a local Coulomb repulsion: 

k,(T 

-2JS^4 - 2J'(5o+,s„- + SoS+) . (4) 

We solve this model using the numerical renormalization 
group method [2^ , which allows calculation of dynamical 
quantities on the real energy axis at both zero and finite 
temperature The impurity self-energy E(a;) is 

calculated using the method described in Ref. [s^l- For 
comparison with ED, it is then evaluated on the imagi- 
nary axis z ~ iLU hy analytic continuation 

= / duj' ^^ . (5) 

The calculations use a logarithmic discretization of the 
conduction band Sk^ — > ±I?2A ~ with A = 1.5 and we 
retain of order 600 states per NRG iteration. Details of 
the calculation of spectra and other dynamical quantities 
can be found in Refs. [Ull^. 

The effective model specified in Eq. (O allows the non- 
Fermi-liquid physics of the intermediate phase of the two- 
band model to be understood qualitatively [Ul . The case 
of Ising exchange (J' = 0) describes interacting conduc- 
tion electrons scattering from a disordered static config- 
uration of local spins with Sf = ±1/2. The disorder 
potential is proportional to the Ising exchange J. The 
conduction electrons therefore have a finite lifetime at 
the Fermi level even at T = 0, as is characteristic of 
a disordered metal. Switching on the spin-flip part of 
the Hund's exchange gives the disordered spin configura- 
tion some dynamics at finite temperature, but as long as 



J' < J, we expect from the ferromagnetic Kondo model 
that the spin-flip part of the Hund's exchange will renor- 
malize to zero at low temperature with the Ising part of 
the Hund's exchange remaining finite. At T = 0, the 
system will again be disordered with a disorder poten- 
tial given by the finite Ising part of the Hund's exchange. 
This results in a finite lifetime for the conduction elec- 
trons at T = 0. 

The situation changes for isotropic Hund's exchange 
(J' = J). In this case, the ferromagnetic exchange in the 
effective irnpurity model 10} is known to be marginally 
irrelevant [3l|,|32|- Thus, both Ising and spin-flip parts 
of the Hund's exchange renormalize to zero and there 
will be no disorder scattering at T = 0, giving rise to 
an infinite lifetime for the conduction electrons at T = 0. 
The vanishing of the self-energy a.tuj — T — also implies 
that the narrow-band spectral function satisfies a pinning 
condition at w = T = 0. In Sections 5 and 6, these 
qualitative considerations will be seen to be in accord 
with the numerical results from ED for the full two-band 
model and NRG results for the effective one-band model. 



4 METALLIC PHASE 

We begin our discussion of the multi-band ED results 
with the region in which both subbands are metallic, just 
below the Mott transition associated with the narrow 
band. To illustrate the applicability of the ED approach 
this region is very useful since it involves important re- 
distribution of spectral weight between low and high fre- 
quencies that must be captured properly in order to de- 
scribe the Mott transition. 

Fig.0|shows a comparison of the subband self-energies 
for Us — 8 and = 6 in the region close to the lower 
critical Coulomb energy for T — 2.5 meV. For Us ~ 8, 
the Yiiiiujn) vary roughly linearly at low frequencies, in- 
dicating that both bands are metallic. Of course, the 
narrow band is more strongly correlated than the wide 
band. The spacing between the lowest excited state and 
the ground state energy in these calculations is typically 
1 — 2 meV, i.e., evaluation of the self-energies in this low 
temperature range is indeed meaningful. 

The results for = 6 are similar, except that the 
self-energy of the narrow band at U — 2.1 eV is in- 
versely proportional to indicating that the Mott 
transition in this case is located between U = 2.0 eV 
and U — 2.1 eV. In addition, at the lowest Matsubara 
frequencies the self-energy of the wide band at 2.1 eV 
exhibits deviations from linear icOn variation. These de- 
viations are related to finite-size effects and will be ana- 
lyzed in more detail in the next section when we discuss 
the intermediate orbital-selective phase. Just below the 
transition for Us — 6, at U = 2 eV, very small deviations 
from approximately linear icun variation can also be seen 
in both Si(?a;„). 
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FIG. 4: Imaginary part of subband self-energies in metallic 
region for J = J' = i7/4 at T = 2.5 meV. Solid (red) curves: 
narrow band, dashed (blue) curves: wide band. From top: 
U = 1.8, 1.9, 2.0, 2.1 eV. (a) n, = 8; (b) n, = 6. 



FIG. 5: Imaginary part of subband self-energies in in- 
termediate phase, for U = 2.4 eV and J — J' — 11/4. 
T = 2.5, 5, 10, 20, 30 meV (from bottom), (a) ris = 8; 
(b) Us = 6; 



It is remarkable that the ED resuhs for Us — 8 capture 
the correlation effects in the common metallic phase of 
both subbands very well until close to Uci and down to 
rather low temperatures. This result is not at all obvious 
since so close to the Mott transition a large fraction of the 
spectral weight of the narrow band is transfered from the 
Fermi level to the Hubbard bands. Even the results for 
Us = 6 are qualitatively correct. The main effect of the 
cluster size = 6 is the underestimate of Uri by about 
0.1 eV. According to the single-band results [16j inclusion 
of the zero energy bath levels in the Ug — 8 cluster should 
provide the most important part of the shift towards the 
correct Ud- In addition these extra bath levels yield an 
excellent representation of the metallic properties of both 
subbands right up to the critical Coulomb energy. 

Of course, slight finite-size effects should be manifest 
also in the low-temperature ED results for Us — 8. The 
precise form of 'Ei{iujn) at very low frequencies and tem- 
peratures, in particular, the range of true Fermi-liquid 
behavior in the immediate vicinity of Ud, can only be 
investigated by more accurate methods, such as a full 
two-band extension of the NRG approach which is appli- 
cable at zero and finite T [ssj . 

We emphasize that this region of the T/U phase di- 
agram for isotropic Hund's exchange is not yet accessi- 



ble using QMC calculations. It would therefore be of 
great interest to extend the present results to Ug = 9 and 
ns = 12 in order to explore the strongly correlated metal- 
lic phase of realistic three-band materials which have so 
far been investigated only for Ising-like exchange inter- 
actions. A more detailed comparison of the role of Hund 
vs. Ising exchange treatments in the metallic phase of 
the present two-band model will be given elsewhere [36| . 

5 INTERMEDIATE PHASE: ISOTROPIC 
HUND'S EXCHANGE 

As shown first by Koga et al. (1J| within ED/DMFT 
calculations at T = 0, the two-band Hubbard model 
with full Hund's exchange interaction exhibits successive, 
orbital-selective Mott transitions. In Ref. |3| we proved 
that analogous ED calculations at finite temperature are 
consistent with these findings, revealing sequential first- 
order transitions associated with the two subbands. In 
this section we discuss in more detail the low temper- 
ature properties of the intermediate phase between the 
fully metallic and insulating phases. 

Fig. 13 shows the subband self-energies for a Coulomb 
energy U = 2.4 eV, i.e., between the lower and upper 
Mott transitions. The results for Us = 8 demonstrate 



7 




s 

W 





-0.2 
-0.4 
-0.6 
-0.8 
-1 



1 

t 


1 1 1 
(a) 




NRG .y.,js^*** 


\ 








1 


ED 

1 1 1 



FIG. 6: Imaginary part of subband Green's functions in in- 
termediate phase, for U = 2.4 eV, J = J' = 11/4, Us — 8. 
T = 2.5, 5, 10, 20, 30 meV. The dot indicates the pinning 
condition for the wide band. 



that the narrow band is insulating, i.e., its self-energy is 
inversely proportional to itOn at low frequencies. T,i(iujn) 
is also nearly independent of temperature in the range 
T = 2.5 ... 30 meV, suggesting that the Mott gap is con- 
siderably larger. The self-energy of the wide band, how- 
ever, extrapolates to zero at low frequencies, indicating 
metallic behavior. These results suggest that the pin- 
ning condition iV2(0) = ^2(0) is satisfied, i.e., the spec- 
tral density at Ep of the interacting system coincides 
with that of the non-interacting system. Because of the 
rapid increase of lniT,2{iuJn) at small ujn, the spectral 
distribution of the wide band in the intermediate phase 
consists of a very narrow peak at Ep and large upper and 
lower Hubbard bands. Nevertheless, as will be discussed 
further below, the comparison with the NRG results in- 
dicates that this band does not satisfy true Fermi-liquid 
behavior. 

The results for Ug ~ 6 shown in Fig. [S^b) agree qual- 
itatively with those for = 8. The self-energies of the 
narrow band are nearly identical, with only a slightly 
larger dependence on temperature in the case rij = 6. 
The main difference is that E2(iw„) at low T and small 
iun exhibits larger deviations as a result of finite-size ef- 
fects. This is plausible since, as discussed above, the wide 
band is metallic, but highly correlated. The pronounced 
three-peak structure of its spectral function is therefore 
not represented accurately by only two non-zero bath lev- 
els per impurity orbital in the ED calculation for ti^ = 6. 
The extra bath level with zero energy included for Us = 8 
therefore plays a crucial role for a proper description of 
the metallic behavior of the wide band for full Hund's 
coupling. 

Fig.|Blshows the subband Greens functions Gi{iuJn) in 
the intermediate phase at various temperatures. ImGi 
is linear in con at low frequencies, as expected for the 
insulating subband. Although the wide band exhibits 
clear signs of finite-size effects at small w„ it is neverthe- 
less seen to approximately satisfy the pinning condition 
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FIG. 7: Imaginary part of self-energy E2 of wide band 
for U = 2.4 eV and J = J' = (7/4. (a) Sohd (red) 
and dashed (blue) curves: ED results for Us = S at T — 
2.5, 5, 10, 20, 30 meV (from bottom); green curves: NRG 
results for T = 3, 10, 15, 34 meV (from top), (b) Same 
quantities on expanded low-frequency scale. Black curve (-I-): 
self-energy Ei of narrow band at 2.0 eV. 



ImG2(«(^n) ^ — 7rA'^2(0) = —1 in the limit LUn — *■ 0. 

For computational reasons the role of finite-size effects 
in the present finite-temperature ED approach based on 
full matrix diagonalization can at present not be checked 
by going beyond the cluster size = 8. To analyze 
the low-frequency properties of the wide band in the in- 
termediate phase more closely, we have carried out NRG 
calculations for the effective one-band model discussed in 
Section 3. Although this model is more appropriate for 
the anisotropic exchange treatment discussed in the next 
section, we use it here (with some caution) since it is so 
far the only scheme capable of providing a guideline for 
the two-band Hubbard model at low finite temperatures. 

In Fig. Efa) the ED results for lmT,2{iuJn) are com- 
pared with the corresponding self-energy derived within 
the NRG. The overall frequency variation of E2 is seen to 
be remarkably similar for both methods, implying similar 
spectral distributions. The most noticeable difference is 
the larger dependence on temperature in the case of the 
ED results. According to Fig. the increase in cluster 
size from Us — 6 to ng = 8 diminishes the variation of E2 
with temperature and makes the minimum near a;„ « 0.9 
less deep. This trend indicates that larger cluster sizes 



8 



0.3 



3 



' 1 ' 


1 ' 1 ' 




j\ NRG 




/ \ ^'"^ 




I.I. 



-4 -2 2 4 

CO (eV) 

FIG. 8: Spectral distribution of wide band in intermediate 
phase, calculated within NRG approach for U = 2.4 eV, J = 
J' = U/4:. Green curve: T = 3 meV; red curve: T — 30 meV; 
dashed blue curve: bare density of states. 



would bring the ED results into better agreement with 
the NRG data. We point out, however, that part of the 
difference should be caused by the approximate nature 
of the NRG model in the case J' = J. As shown in the 
following section, for J' = the ED and NRG results for 
S2 at T = 20 ... 30 meV are in perfect agreement. 

The low- frequency region of E2(i(j;n) is plotted in 
greater detail in Fig. Cfb). In the ED calculation, the 
spacing between the lowest excited level and the ground 
state energy is typically 10 meV. Thus, temperatures be- 
low this range will give rise to increasing finite-size effects, 
precluding any meaningful analysis at low frequencies. 
Despite this limitation, the ED data follow the frequency 
dependence of the NRG results rather well. 

In order to provide some insight into the variation of 
Ti2 at finite aj„, we compare it in Fig. [Tfb) also with the 
self-energy of the narrow band Si at f/ = 2.0 eV, i.e., just 
below the Mott transition at Uci « 2.1 eV. According to 
Fig. 1(a), we find: Zi{U = 2.0 eV) « Z2{U = 2.4 eV) « 
0.11. Nevertheless, whereas Ei at [/ = 2.0 eV exhibits 
the typical behavior expected for a strongly correlated 
Fermi-liquid phase, S2 evidently has much larger nonlin- 
ear corrections. A similar systematic difference is found 
between E2 in the intermediate phase and the behavior 
the self-energy in a single-band model for U just below 
the Mott transition. 

The NRG results shown in Fig. indicate that 
'E2(ioJn) ^ in the limit Wn —>^ and for low T. Thus, 
particles at Ep have infinite lifetime. In addition, anal- 
ogous NRG calculations at T = yield 1351 non- Fermi- 
liquid behavior, in agreement with Ref. |2l| . Thus, at 
small finite iujn Im E2 does not vary linearly. The quali- 
tative agreement between the finite-T NRG and ED re- 
sults seen in Fig.[7|shows that the ED data are consistent 
with non-Fermi-liquid behavior. 

Fig. |S1 shows the spectral function of the wide band 
in the intermediate phase, calculated within the NRG. 



Spectra of this kind were also obtained by Arita and Held 
[g for U near Uci within QMC/DMFT calculations for 
the same two-band model for T — and J' ~ J ~ U/A. 
As mentioned above, the distribution exhibits a very nar- 
row peak which satisfies the pinning condition at Ep for 
temperatures up to about 10 meV. Because of the sharp 
central peak of A2{uj), the real and imaginary parts of the 
self-energy of the wide band Tt2{(jj) exhibit structure on 
a similar scale. Nevertheless, the variation of ImE2(ia;„) 
obtained in the NRG is very smooth. Thus, the weak 
shoulder seen in the ED results for S2 near u;„ « 0.06 
in Figs.|3ta),IHa) and|7|must be attributed to finite-size 
effects. 

The ED and NRG results discussed above demonstrate 
that, in the presence of full Hund's exchange, the wide 
band above the main Mott transition at Uci remains 
metallic, without satisfying Fermi-liquid criteria. In the 
following section it will be shown that the absence of 
spin-flip and pair-exchange terms enhances this trend to- 
wards non-Fermi-liquid behavior, so that even particles 
at Ep acquire a finite lifetime. Thus, in both cases, the 
intermediate phase is bad-metallic. 

We point out that the differences between the self- 
energies for Us — 8 and = 6 appear mainly at rather 
low temperatures, below about T = 20 meV. Since 
QMC calculations also become difficult to converge in 
this range, three-band ED calculations with two bath lev- 
els per impurity orbital, i.e., cluster size rig = 9, should 
be competitive with analogous QMC calculations - with 
the important additional benefit, that the ED approach 
can handle full Hund's coupling. 

6 INTERMEDIATE PHASE: ISING EXCHANGE 

As shown in Ref. ,13], in the absence of spin-flip and 
pair-exchange contributions, the first-order Mott tran- 
sition at Uci ~ 2.1 eV affects the different subbands of 
the Hubbard model in fundamentally different ways: The 
narrow band undergoes a complete metal insulator tran- 
sition, but the wide band changes from a normal metal to 
a bad metal in the sense that its self-energy exhibits pro- 
gressive deviations from Fermi-liquid linear cOn variation 
at low frequencies. Thus, Im T,2{ii-Un) — > c{U) for ujn — >■ 0, 
where c{U) w ^ —00 for U — Uci Uc2- Spec- 
tral functions obtained via the maximum entropy method 
showed that this breakdown of Fermi-liquid behavior in 
the wide band in the intermediate phase leads to a nar- 
row pseudogap near Ep and to a violation of the pinning 
condition, i.e., ^2(0) < iV2(0). The spectral function of 
the wide band then acquires a characteristic four-peak 
structure, with two maxima flanking the pseudogap, in 
addition to the Hubbard bands at higher energies. With 
increasing U , the pseudogap becomes deeper and wider, 
until it turns into a true insulating gap at Uc2 ~ 2.7 eV. 
As also shown in Ref. |l3l |. the upper bad- metal to insu- 
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FIG. 9: Imaginary part of subband self-energies T.i{iuj„) 
in intermediate phase for Ua = 8, U = 2.4 eV, J = (7/4, 
J' = 0; Solid (red) and dashed (blue) curves: ED results for 
T = 2.5, 5, 10, 20, 30 meV (from bottom), (a) aj„Ei(iaj„); 
(b) T,2{icL}n); green curves: NRG results as in (c). (c) Low- 
frequency behavior of self-energy of wide band: red and blue 
curves: ED results as in (b); green curves: NRG results for 
T = 2, 4.4, 10, 20, 30 meV (from bottom). For clarity, suc- 
cessive curves in (c) are displaced vertically by 1. 



lator transition is not a first-order transition, in contrast 
to the main transition at Ud- 

In this section we discuss the temperature dependence 
of the subband self-energies in the intermediate phase 
and the finite-size effects associated with the limited 
number of bath levels included in our ED/DMFT ap- 
proach. 

Fig. |51 shows the subband self-energies for Ug = 8 and 
U = 2.4 eV, assuming again J = f//4, but J' = 0. The 



FIG. 10: Imaginary part of subband self-energies T,i{iLu„) in 
intermediate phase for Us = 6, U = 2.2 eV and J = (7/4, 
J' = 0. T = 2.5, 5, 10, 20, 30 meV (from bottom). 



narrow band is fully insulating, so that uJn IniSi(«ti-'n) — 
const, nearly independently of temperature, similarly to 
the case J' ^ J plotted in Fig. |S1 The self-energy of the 
wide band, however, differs qualitatively from the behav- 
ior found for isotropic Hund's coupling. Rather than van- 
ishing in the limit w„ — > 0, Iml]2(iwn) now approaches 
a sharp minimum, which gets more pronounced towards 
low temperatures. Clearly, since w„ImE2(i^ri) —> at 
small LUn, this band is not yet insulating. Thus, instead 
of satisfying the pinning condition at Ep, the spectral 
function of the wide band exhibits a dip or pseudogap at 
the Fermi level (see below). 

Fig. Efb) also shows the NRG results for 'S2{'ii^n)- 
They are seen to be in very good agreement with the 
ED data, except at low frequencies for T < 10 meV. The 
fact that for T = 20 ... 30 meV there is now much bet- 
ter coincidence between the ED and NRG results than 
in Fig. [7| for ,]' = ,] indicates that, as discussed in Sec- 
tion 3, the effective one-band model employed for the 
NRG is more appropriate in the case J' — Q. The low- 
frequency behavior of Im S2(iw„) is shown in more detail 
in Fig. ^c). Both the ED and NRG results reveal that 
the sharp minimum of E2 at small ojn gets progressively 
deeper at low temperatures, suggesting that the pseudo- 
gap in the spectral function becomes accordingly deeper. 

Analogous results for = 6 are given in Fig. 1101 Since 
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FIG. 11: Imaginary part of subband Green's functions in 
intermediate phase, for U = 2.4 eV, J — U/4, J' — 0, Us — 8. 
T = 2.5, 5, 10, 20, 30 meV (from top). 



FIG. 12: Spectral distribution of wide band in intermediate 
phase, calculated within NRG approach for U = 2.4 eV, J = 
U/4:, J' = 0. Green curve: T = 3 meV; red curve: T = 
34 meV; dashed blue curve: bare density of states. Black 
dotted curve: spectrum for J' = J, T = 3 meV from Fig. |S| 



according to Ref. |16| in this case the wide band becomes 
insulating near U ~ 2A eV, we choose U = 2.2 eV to il- 
lustrate the non- Fermi-liquid behavior of this band in the 
intermediate phase. As for Us — 8, the narrow band is in- 
sulating, i.e., uJnlm'E,i{iujn) const, at low frequencies. 
The frequency variation of 'E2{i^n) is seen to be more 
strongly affected by finite-size effects. Nevertheless, as in 
the case Ug — 8, a sharp minimum is found for a;„ 0. 

Fig. II II shows the subband Green's functions Gi{iujn) 
for anisotropic Hund's exchange at various temperatures. 
As for the case J' = J shown in Fig.El ImGi is linear in 
ujn at low frequencies, since the narrow band is insulating. 
On the other hand, because of the more severe breakdown 
of Fermi-liquid behavior in the wide band for J' = 0, with 
IniE2(i^n) 7^ in the low-frequency limit, the Green's 
function no longer satisfies the pinning condition. Thus, 
IniG2(it^„) -c{T) with c{T) < 7riV2(0) = 1, and 
c(r) — * for decreasing temperature. This behavior im- 
plies that the spectral function of the wide band exhibits 
a pseudogap at Ep which becomes progressively deeper 
towards low T. 

This picture is fully confirmed by the NRG results 
shown in Fig. 1121 In contrast to the narrow peak at Ep 
in the three-peak structure seen in Fig. |S1 for isotropic 
Hund's coupling, the spectra for J' = show a pseu- 
dogap which becomes deeper as T decreases. As a re- 
sult, the spectral distribution now exhibits a character- 
istic four-peak structure, with two maxima limiting the 
pseudogap and two shoulders associated with the Hub- 
bard peaks. Spectra of this kind were first observed in 
the QMC/DMFT results at T = 31 meV reported in 
Ref. [l3j (see also next section) . The comparison with the 
spectra for J' = J reveals nearly identical excitations for 
\lo\ > J. Thus, as expected, the diff'erent treatments of 
exchange interactions affect primarily the low-frequency 
excitations in the metallic wide band. 



7 COMPARISON WITH PREVIOUS 
QMC/DMFT RESULTS 

In this and the following sections we compare the ED 
and NRG results for J' = with available QMC/DMFT 
data in order to illustrate the consistency between these 
impurity treatments and to explore further the role of 
finite-size effects in the ED approach. 

As shown in Ref. i^i the absence of spin-flip 

and pair-exchange terms, purely metallic and insulating 
phases exist for U <Uci~2.l eV and U >Uc2~ 2.7 eV, 
respectively. Fig. 1131 shows a comparison of subband 
self-energies at three representative Coulomb energies 
within the intermediate 'bad-metal' non-Fermi-liquid re- 
gion. The ED results for = 8 are seen to be in excel- 
lent agreement with the NRG self-energies, confirming 
the validity of the effective one-band model in the inter- 
mediate phase for J' — 0. Moreover, both schemes agree 
very well with the QMC/DMFT self-energies reported in 
Ref. for T = 31 meV. Minor differences between the 
QMC and ED/NRG results are found only in the steep- 
est parts of Si close to the first Matsubara frequency. (In 
contrast to the QMC calculations which are carried out 
at discrete Matsubara frequencies, the NRG self-energy 
is available continuously as a function of frequency. The 
ED self-energies could in principle also be obtained at ar- 
bitrary iLU but were calculated here at iuj„ . This explains 
the slightly different form of the curves plotted in panel 
(c). At iujn the ED and NRG data nearly coincide.) 

The overall trend obtained previously within the 
QMC/DMFT is fully confirmed by the new ED and NRG 
calculations: The narrow band is insulating throughout 
this range of Coulomb energies, whereas the wide band 
changes gradually from metallic to insulating via pro- 
gressive non- Fermi-liquid behavior. Thus, for small u!„, 
Iml]2(iaj„) — > c{U), where c{U) -oo at U ^ 2.7 eV. 
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FIG. 13: Comparison of subband self-energies Ei(jcj„) in in- 
termediate non-Fermi-liquid phase for U = 2.1, 2.4, 2.7 eV, 
J = J7/4, J' = 0, r = 31 meV, calculated via three different 
impurity solvers: (a) ED results for = 8: red solid curves: 
narrow band; blue dashed curves: wide band; (b) QMC re- 
sults from Fig. 10 of Ref. p^ : red solid curves: narrow band; 
blue dashed curves: wide band; (c) comparison of ED (x) and 
NRG self-energies of wide band. 



Evidently, despite their intrinsic numerical uncertain- 
ties, the three complementary impurity solvers: ED, 
NRG and QMC provide perfectly consistent descriptions 
of the electronic properties of the orbital-selective phase 
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FIG. 14: Subband self-energies Ei(ia;„) in intermediate non- 
Fermi-liquid phase for T = 31 meV, U = 2.1, 2.4, 2.7 eV, 
J = U /4, J' = 0. Solid (red) curves: Us — 8; dashed (blue) 
curves: Us = 6. (a) narrow subband: aj„Ei(iaj„); (b) wide 
subband: E2(iti;„). 



for Ising-like exchange. 

Since the temperature T = 31 meV in these results 
is relatively high, finite-size effects tend to be less pro- 
nounced than in the cases discussed in the previous sec- 
tions at lower T. This is illustrated in Fig. 1141 where we 
compare the above ED results for Ug = 8 with those for 
Us = 6. The self-energies of the narrow band are nearly 
identical and satisfy ij;„Ei(ia;„) — > const, at low frequen- 
cies. This is plausible since the additional zero-energy 
bath level for — 8 carries very little weight in the 
insulating state. Because of breakdown of Fermi-liquid 
behavior, the self-energy of the wide band approaches a 
finite value in the limit — > 0. The differences between 
the results for Ug = 6 and Hs = 8 are very small as long as 
this band is either nearly metallic, like at U = 2.1 eV, or 
nearly insulating, like at [/ — 2.7 eV. Slightly larger dif- 
ferences are found only in the middle of the bad-metallic 
region near U = 2A eV. As discussed in the preceding 
section, the spectral function then has a more compli- 
cated four-peak structure as a result of the narrow pseu- 
dogap at Ep, which cannot be adequately represented 
using only two bath levels. This region could possibly be 
even more accurately described by using four bath lev- 
els for the wide band and two levels for the insulating 
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FIG. 15: Spectral distributions of subbands in intermedi- 
ate phase, calculated within QMC/DMFT for U = 2.4 eV, 
J = C//4, J' = 0, T = 31 meV. Red curve: metallic wide 
band; magenta curve: insulating narrow band; blue and black 
curves: bare densities of states. From Fig. 11(b) of Ref. [l^. 

narrow band, maintaining the total Ug — 8. 

We close this section by showing in Fig.^|the spectral 
distributions of both subbands, as calculated within the 
QMC/DMFT and the maximum entropy method ^13,J . 
The spectrum of the bad-metallic wide band can be com- 
pared with the corresponding NRG spectra plotted in 
Fig. El Both distributions exhibit the marked four-peak 
structure induced by the pseudogap and the Hubbard 
bands. The low-frequency region is seen to be in ex- 
cellent agreement. Both methods coincide in that, at 
T « 30 meV, the interacting density of states at Ep 
of the wide band is only about one third of the nonin- 
teracting one. The slightly larger differences at higher 
energies, in particular, the position and relative weight 
of the Hubbard bands, may be caused by the choice of 
fitting parameters in the maximum entropy procedure, 
and by the less accurate nature of the NRG approach at 
high energies. 



8 COMPARISON WITH QMC/DMFT RESULTS 
FOR W2 = lOWi 

The breakdown of Fermi- liquid behavior in the orbital- 
selective phase with Ising-like exchange was recently also 
studied by Biermann et al. for the same two-band 
Hubbard model as above, except for Wi = 0.2 eV and 
W2 = 2 eV. The subbands were also assumed to be non- 
hybridizing and half-filled. QMC/DMFT calculations 
were carried out for T = 1/120 . . . 1/40 eV, U = 0.8 eV, 
J = and J = U/4:, with J' — 0, i.e., in the absence of 
spin-flip and pair exchange terms. For J = as well as 
J — C//4, the self-energy of the narrow band was found 
to diverge at low frequencies, demonstrating that this 
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FIG. 16: Zi{U) as a function of U for Wi = 0.2 eV, W2 = 
2.0 eV, T = 1/120 eV, calculated within ED/DMFT, Us = 6, 
for J = J' = 0, J = J' = f//4, and J = U/4:, J' = 0. Red and 
blue curves: wide band, black curves: narrow band. The Zi 
for J — J' — U/4: and J = U/4, J' = are indistinguishable. 
The dots mark the three cases shown in Fig. 1171 

band is insulating since Wi ^ U. The wide band is 
in the intermediate orbital-selective region between the 
purely metallic and insulating phases since U ^ W2 ■ The 
self-energy of this band showed a striking variation with 
the magnitude of J: For J = 0, Tt2{iuJn) at low 

frequencies, as expected for Fermi-liquid behavior. For 
J = U/4:, however, ImE2(iw„) const, w —0.09 nearly 
independently of temperature for T — 1/1 20. ..1/40 eV, 
indicating breakdown of Fermi-liquid properties. Accord- 
ingly, the pinning condition iV2 (0) — A2 (0) was found to 
be satisfied for J = 0, but not for J = U/4. These results 
are fully consistent with the trend discussed in Ref. 0| 
for 14^1 2 eV and W2 = 4 eV. 

To check the accuracy of our finite-T ED approach 
we have applied it to the case investigated in Ref. 21]. 
To provide a picture of the Mott transitions in this 
two-band system we show first in Fig. 1161 the variation 
of Zi{U) for three different treatments of Hund's ex- 
change. In all cases the narrow band becomes insulating 
at about Uci = 0.3 . . . 0.4 eV. However, the range and 
nature of the orbital-selective phase of the wide band, 
Uci < U < Uc2, depend sensitively on the magnitudes of 
J /U and J' /U . For J = J' = the upper transition oc- 
curs for Uc2 slightly larger than W2, with a pronounced 
hysteresis loop indicative of first-order behavior. A very 
weak hysteresis is found also for J = J' = U/4, with 
Uc2 ~ 1-5 eV, similar to the one in Fig. linear J7 « 3 eV. 
Finally, for J = U/4 and J' = 0, C/c2 w 1.2 eV with- 
out evidence of first-order behavior. (The latter result 
is consistent with Uc2 ~ 2.4 eV in Fig. 2 of Ref. 0.) 
Although Z2{U) « 0.2... 0.8 at U ^ 0.8 eV, i.e., for 
Uci <C i7 <C Uc2, the analysis of the self-energy reveals 
fundamentally different electronic properties of the wide 
band in this region. 

This is illustrated in Fig. 1171 which shows Iml]2(*w„) 
for Us = 8 and Us = 6- The ED results for = 8 
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FIG. 17: (a) Imaginary part of self-energy of wide band for 
J = and J = [7/4 in intermediate phase for U = 0.8 eV, 
J' = 0, calculated within ED/DMFT with =8. Wi = 
0.2 eV, W2 = 2.0 eV. Red solid, blue dashed, and black dotted 
curves: T = 1/120 eV, T = 1/80 eV, and T" = 1/40 eV, 
respectively. The dots denote the results for T = 1/120 eV 
obtained in Ref. f2l|. (b) Same as (a) except for Us = 6. (c) 
Same as (a) except ior J' = J = U /4. 



are in excellent agreement with the QMC data |2l| for 
J = J' = as well as J = C//4, J' = 0, with the excep- 
tion of small finite-size effects at the lowest Matsubara 
frequencies and lowest temperature. (Since W2 = 2 eV, 
T = 1/1 20. ..1/40 eV in these calculations corresponds 
to about T = 16 ... 50 meV in the case W2 = 4 eV con- 
sidered in the preceding sections.) Also, the variation 
with T is in our results slightly more pronounced than in 
the case of the QMC/DMFT. 

For Us = 6 the ED results exhibit larger finite-size ef- 



fects, as shown in Fig. llTr b). The variation with temper- 
ature is in this case also larger than for = 8. Evidently, 
the subtle features of spectral functions in the interme- 
diate phase are not so well represented by including only 
two bath levels. Nevertheless, the qualitative difference 
between Fermi-liquid behavior for J = and the clear 
deviation from this behavior for J = J7/4 are very well 
reproduced by these tis = 6 ED/DMFT calculations. 

For completeness we show in Fig. [TTT c) the — 8 
ED/DMFT results for isotropic Hund's coupling, i.e., 
J' ~ J — U/A. This case is not yet accessible within 
QMC /DMFT because of sign problems at low tempera- 
tures. The inclusion of spin-flip and pair-exchange terms 
is seen to restore the limiting behavior T,2(iujn) —>- for 
iVn 0. Nevertheless, compared to the case J' = J = 
shown in (a), S2 now increases much more rapidly at 
small finite frequencies, similarly to the data shown in 
Figs. [SI and [71 (As argued in Section 5, the shoulder near 
ujn ~ 0.06 is caused by finite-size effects.) Thus, instead 
of quasiparticles with weight Z2 ~ 0.8 for J = J' = 0, 
finite Hund's exchange with J' ^ J = [7/4 supports a 
state of infinite lifetime at Ep , but does not satisfy Fermi- 
liquid criteria at > 0. This picture is consistent with 
the T = ED/DMFT results by Biermann et al. [2l|. 
The spectral function of the wide band therefore should 
exhibit a narrow peak at Ep, similar to the one shown 
in Fig. [HI for U = 2.4 eV. 



9 COMMENTS ON IPT/DMFT RESULTS 

In Ref. ^3 reported DMFT calculations for the 
two-band Hubbard model within QMC and iterated per- 
turbation theory [s^. Both approaches were shown to 
yield a consistent picture of the Mott transition, in the 
sense that both exhibit a single first-order transition 
at which the narrow band becomes insulating and the 
wide band begins to show progressive bad-metallic be- 
havior. This band becomes fully insulating at a higher 
Coulomb energy in a non-first-order manner. As shown 
via ED/DMFT at T = in Ref. |li| and at T > in 
Ref. [1^, the nature of this upper transition depends in 
a critical way on the treatment of onsite exchange in- 
teractions. In the absence of spin-flip and pair-exchange 
terms it is continuous, consistent with the QMC results, 
while for full Hund's coupling it becomes first-order. 

The question then arises whether the IPT approach 
also supports this picture. In fact, as pointed out in 
Ref. [13, it is surprising that the IPT gives only one 
common first-order transition even though spin-flip and 
pair-exchange terms are included. It can easily be shown, 
however, that to second-order in the Coulomb interaction 
the subband self-energies with and without spin-flip and 
pair-exchange are the same, except for a slightly different 
relation between U and J. Thus, J' — J = 0.25 U gives 
the same result as J' = and J « 0.22 U. Corrections to 
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the simple second-order diagram via renormalization of 
subband energies to yield the correct atomic limit [s^ls^ 
are also insensitive to the choice of J' for the present two- 
band model. 

To resolve this puzzle we have repeated the IPT calcu- 
lations reported in Ref. at even lower temperatures 
and found indeed a tiny hysteresis loop also at the up- 
per Mott transition, with a critical temperature of ap- 
proximately ~ 5 meV, i.e., significantly lower than 
the critical temperature of the lower Mott transition, 
Tci w 50 meV. These findings suggest that present formu- 
lations of IPT are too simple to deal with the full com- 
plexity of Hund's exchange. Diagrams beyond second- 
order are required to distinguish more clearly between 
J' — J and J' — treatments. 



10 SUMMARY AND OUTLOOK 

Finite-temperature ED /DMFT studies are carried out 
in order to explore the usefulness of this approach for 
multi-band systems. The important feature here is that 
onsite Coulomb and exchange interactions are fully in- 
cluded. Since in the past ED has been used mainly for 
single-band systems, a particular aim of this work is to 
illustrate the dependence of self-energies on the cluster 
size and to test the range of applicability of this method. 
As an example, we focus on the Hubbard model for two 
bands of different widths and investigate the metal insu- 
lator transition as a function of Hund's exchange. 

The surprising and potentially very useful result of this 
work is that even a cluster size of rig = 6, with only two 
bath levels per impurity orbital, provides a qualitatively 
correct picture in all of the important phases of the T/U 
phase diagram. Moreover, finite-size effects are substan- 
tially reduced for Us — 8, i.e., using one extra bath level 
per impurity orbital. Thus, in the phase just below the 
main first-order Mott transition, both subbands exhibit 
clear metallic properties, albeit with strongly reduced 
quasi-particle weights. The intermediate region between 
the purely metallic and insulating phases depends in a 
subtle manner on the exchange interactions included in 
the ED calculation. For full Hund's coupling coexist- 
ing metallic and insulating subbands are found, where 
the wide band exhibits infinite lifetime at Ep, but non- 
Fermi-liquid behavior at finite frequencies. For Ising-like 
exchange, this breakdown of Fermi-liquid behavior is en- 
hanced, giving finite lifetime even at Ep- 

To explore the influence of finite-size effects in the two- 
band ED/DMFT approach NRG calculations were per- 
formed for an effective one-band model suitable to de- 
scribe the wide band in the intermediate phase. Both for 
full Hund's coupling and Ising-like exchange, very good 
agreement with the ED results is found for temperatures 
in the range T = 20 ... 30 meV. At lower T, larger dif- 
ferences appear on a quantitative level at low frequen- 



cies. Nevertheless, the important qualitative differences 
between various phases, especially the characteristic low- 
frequency variation of the self-energy in the two types 
of non-Fermi-liquid regions for isotropic and anisotropic 
exchange coupling, are fully reproduced by the ED ap- 
proach, both for Us — 6 and = 8. 

To test the accuracy of the two-band ED approach 
we also have applied it to models studied earlier within 
QMC/DMFT, neglecting spin-flip and pair-exchange. 
Nearly quantitative agreement is obtained for a cluster 
size Us = 8, but even the ED results for = 6 are found 
to be qualitatively reliable. 

The present ED approach utilizes full diagonalization 
of the impurity Hamiltonian. Since at low temperatures 
only a limited range of excited states is relevant for the 
local Green's functions and self-energies it should be very 
useful to generalize finite-T Lanczos one-band methods 
[2^ in order to bridge the gap between the present 
work and the true T = limit, and to apply the finite- 
temperature ED/DMFT approach to realistic two-band 
and three-band materials. 
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